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Abstract. For the class of quasi-periodic solutions of the vortex filament equation we study 
connections between the algebro-geometric data used for their explicit construction, and the 
geometry of the evolving curves. We give a complete description of genus one solutions, 
including geometrically interesting special cases such as Euler elastica, constant torsion 
curves, and self-intersecting filaments. We also prove generalizations of these connections 
to higher genus. 

1. Introduction 
In 1972, H. Hasimoto J7j introduced the complex function 

g=|/texp^i J rds*j (1) 

of the curvature k and the torsion r of a space curve, and showed that if the position vector 
7 of the curve evolves according to the vortex filament equation (VFE), 

fry = fry d 2 7 

dt ds <9s 2 ' 1 ' 

then the potential q(s,t) is a solution of the focussing cubic nonlinear Schrodinger equation 
(NLS) 

iq t + q ss + 2\q\ 2 q = 0. (3) 

In these equations, s is the arclength parameter of the curve and the spatial variable for 
the NLS potential, and t is time. In fact, the evolving curve defines q(s,t) uniquely, up to 
multiplication by a unit modulus constant. 

The association *y(s, t) > q{s,t) is often referred to as the Hasimoto map. Since the NLS 
had been shown to be completely integrable by Zakharov and Shabat El- Hasimoto's 
discovery meant that the VFE is also completely integrable, and possessed all the attendant 
structure: infinite sequences of commuting flows and common conservation laws, solution 
by inverse scattering, and special solutions such as solitons and, more generally, finite-gap 
solutions. 

The VFE was proposed a century ago as a simplified model for the self-induced motion 
of vortex lines in an incompressible fluid |3*T] I25j. Hasimoto's discovery has developed into 
one of the richest examples of the connection between curve geometry and integrability. The 
understanding of this connection has progressed in recent years along different directions. In 
the case of the VFE, its bihamiltonian structure and recursion operator have been given a ge- 
ometrical realization j2HlE2] and its relation to the NLS has been shown to possess a natural 
geometric interpretation [2312]. The study of special solutions, including solitons, finite-gap 
solutions and their homoclinic orbits has been undertaken by several researchers, employing 
techniques ranging from perturbation theory fZT\ ITU] to methods of algebraic-geometry and 
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if) ts is equivalent to (jSJ). Expressed in terms of the Pauli matrix 03 = ( \ ^ ) , the matrices 



Backlund transformations [HUE! El, in particular with the aim of understanding whether 
the presence of infinitely many symmetries and integral invariants implies restrictions on the 
topology of the evolving curves |2H| 

In this paper, we focus on connections between geometric information about the evolving 
curve 7, and features of the NLS potential q(s, t). Specifically, we limit ourselves to the class 
VFE solutions that correspond to periodic and quasiperiodic finite-gap solutions of NLS, 
which (as explained in ^2J may be constructed using algebro-geometric data, and it is this 
data which we seek to connect with the geometry (and topology) of 7. Note that this choice 
is particularly significant for closed curves, since the finite-gap solutions are dense in the 
space of all periodic solutions [13], and it is therefore important to study their geometrical 
and topological properties. 

Of course, the class of curves that correspond to finite-gap solutions of the NLS only 
makes sense if the Hasimoto map is invertible. In fact, the inverse map is defined using the 
solutions of the AKNS system for (J3J). As is well-known, the AKNS system consists of a pair 
of first-order linear systems 

l t = vl, ® 

for a vector-valued function t/j, for which the solvability or "zero curvature condition" if) st = 

on the right-hand sides in (@J) are 

^-iA.3 + g ^i(| 9 |*-2AV3 + ( 2iA5 +5< 2 V 

These matrices depend on s and t through the complex-valued function q, and on the spectral 
parameter X. (In equations CQ),© and (0J), and elsewhere in this paper, our conventions have 
been chosen so as to be consistent with the important monograph by Belokolos, Bobenko et 
al. P] . ) Because the first equation of (@J) can be written as Lij) = Xt^ for a first-order matrix 
differential operator L, i/> is known as an eigenfunction associated with q. 

If one can compute a common solution of (@J) for a given potential q(s,t), then the curve 
that corresponds to q under the Hasimoto map is obtained by means of the following recon- 
struction procedure: 

Proposition 1.1 (after A. Sym [H], K. Pohlmeyer [2H]). Let q(s,t) satisfy ©, let V(s,t; X) 
be a matrix whose columns are linearly independent solutions of (@]) ; such that \P(0, 0; A) is 
a fixed matrix in SU{2), and let 

1 d^ 

r(s,t) = V- 1 — . (5) 

dA A=0 

Then T satisfies (0), and the curvature and torsion of T satisfy (JT]) ; assuming we identify 
5u(2) with M. 3 via a linear isomorphism F , specified by 

F: (q 1 fj ~ ei , F: Jj ^ e 2 , F : (° Jj - e 3 . 

(Of course, the reason this construction makes sense is that U, V take value in su(2) when 
A G M.) Thus, if we begin with a solution of 0, construct the potential q(s,t), and solve 
the linear system, then the curve given by © will differ from what we started with by an 
isometry of M 3 . 
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Formula (0) has become known as the Sym-Pohlmeyer reconstruction formula. More 
generally, if in (j5J) we evaluate at an arbitrary real value of A, i.e., 



dA 



Ao e 



(6) 



A=A 



then T still satisfies (j2J), but with a potential that differs from what we started with by 
the transformation q(s,t) i— > e^ as ~ a fS) q(s — 2at,t), a G M, which preserves solutions of (jSJ). 
We will have recourse to the generalized Sym-Pohlmeyer formula © in order to associate a 
closed curve to a periodic NLS potential. 

By differentiating (0) with respect to s, we see that the unit tangent vector T along 7 
and a pair of unit normal vectors U\ , U2 are given by 







Ui = y 



1 
-1 



u 2 = 



with the aforementioned identification of su(2) and M 3 understood. These vectors do not 
constitute a Frenet frame for 7, but rather satisfy the generalized natural frame equations: 



d 



where ki, k 2 are the natural curvatures, which in this case coincide with the real and imagi- 
nary parts of 2q, respectively. (If A = 0, then (T,Ui,U2) comprise a natural or relatively 
parallel adapted frame [2 j along 7.) The interpretation of the Hasimoto map in terms of 
natural curvatures, and the connection between the eigenfunction matrix \1/ and the natural 
frame, will be used when we discuss closure conditions in 
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Organization of the paper. In £0 we introduce the notion of finite-gap solutions, and give 
formulas for the eigenfunctions, potential and curve in terms of the algebro-geometric data — 
specifically a hyperelliptic Riemann surface S and a positive divisor 22 on it, satisfying certain 
reality and genericity conditions. We also use these formulas to give conditions under which 
the curve is smoothly closed, and conditions under which it periodically intersects itself. 
(The closure conditions duplicate results of Grinevich and Schmidt [TSj, but here we derive 
them directly from the curve formulas.) 

In §21 the ideas of the previous section are worked out in detail for the genus one case. 
In particular, we show that the resulting curves are centerlines for Kirchhoff elastic rods. 
We study the geometric features of curves associated with special configurations of the 
algebro-geometric data: we show that the elastic rod centerline becomes an Eulerian elastic 
curve precisely when the zeros of a certain meromorphic differential involved in the finite-gap 
construction coincide, and that the centerlines have constant torsion exactly when the branch 
points of the elliptic curve have an extra axis of symmetry. (They are always symmetric under 
complex conjugation.) 

In £01 we discuss some generalizations of this phenomenon — symmetric branch points 
giving rise to special geometric features of the curve — in higher genus. In an appendix §3] we 
relate our eigenfunction formulas to those given in pQ , and give a straightforward derivation 
of the reality conditions, as well as some new formulas for calculating the algebro-geometric 
data. 
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Forthcoming papers in this series. One of our main interests is the extent to which topological 
information (in particular, knot type) can be "read off" from the algebro-geometric data that 
gives rise to a closed filament. However, it is difficult to obtain examples of closed finite-gap 
filaments to base conjectures on in the first place. As will become apparent in the 
conditions under which the filament is smoothly closed are difficult to compute in general. 

One alternative to tackling the closure conditions head-on is to begin with a closed curve 
and deform its spectrum (in particular, the branch points of the hyperelliptic curve) in a way 
that preserves closure. Deformations that preserve the periodicity of q have been studied 
by Grinevich and Schmidt [T^]. In our next paper, we will examine a special adaptation 
of these deformations that also preserves closure. We will show that, by starting with a 
multiply-covered circle (corresponding to a simple plane wave solution of NLS) and applying 
successive deformations, we obtain a labelling scheme that matches the deformation data 
with the knot type of the resulting filament, assuming the sizes of the deformations are 
sufficiently small. 

2. Finite-gap Solutions 

Soliton equations, when considered on periodic domains, admit classes of special solutions 
which are the analogues of solitary waves for rapidly decreasing initial data on an infinite 
domain. An N-phase, or finite-gap, quasiperiodic potential is an NLS solution of the form 

q(s, t) = q(0 1 , ...,9 N ), 0i(s, t) = kiS + W{t, i = 1, . . . , N, 

such that q is periodic in each phase: q(0i, . . . ,0{ + 2tc, . . . , On) = q(0i, ■ ■ ■ , 0%, ■ ■ ■ , On)- The 
vector k = (ki, . . . , /cat) is called the vector of spatial frequencies, and w = (wx, . . . ,Wn) the 
vector of time frequencies. 

Finite-gap solutions of soliton equations such as KdV, NLS, and general AKNS systems 
were first constructed in the mid-1970's (see the Introduction in pQ for references). In 
particular, finite-gap solutions to the nonlinear Schrodinger equation were first constructed 
by Its ^Hl and Its and Kotlyarov Here, we use Krichever's construction [21] based 

on the Baker- Akhiezer eigenfunction, and the Sym-Pohlmeyer reconstruction formula (j3J) to 
derive the explicit formulas for finite-gap solutions of the VFE (J2J). We mostly follow the 
notation and the discussion of reality conditions in pQ. 

A finite-gap solution is defined by a set of data on a smooth hyperelliptic Riemann surface 
E of genus g, with distinct branch points {Aj}^ 2 , described by the equation 

29+2 

i=l 

Let it : E — ► C U {oo} be the standard hyperelliptic projection defined by it(P) = A, with 
P = (A, fi) G E, and let oo + + oo_ be the pole divisor of n (labelled so that yU" 1 A 9+1 — > ±1 as 
A — ► oo ± ). Let V = P\ + - ■ - + P g+ i be a non-special divisor of degree g + 1 and not containing 
oo±. Recall that a positive divisor T> of degree at least g is called non-special when the linear 
space of meromorphic functions with poles in T> has dimension 1{T>) = deg(X>) — g + 1. The 
special divisors, which are those such that l(T>) > deg(P) — g + 1, are the critical points of 
the Abel map. 

Krichever's main idea is to construct a function ip(P) on E which is uniquely defined by 
a given set of singularities and by prescribed asymptotic behavior near oo ± . 

Definition 2.1. A Baker- Akhiezer function associated to (E, "D, oo±) is defined by the fol- 
lowing properties: 
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• ij; is meromorphic on E\oo± and has pole divisor contained in T>, 

• if) has essential singularities at oo± that locally are of the form 

^ „ e ±Kx s+ 2xh) [ c(S)t) + C ,(A- 1 )] 

for P — > oo± and A = ir(P), where c(s,t) is a function of arbitrary complex 
parameters s and t. 

The singular structure of ip and its normalization at the essential singularities define it 
uniquely as described in the following 

Proposition 2.2 (Krichever Suppose that the following technical condition holds: 

Genericity Condition: The non-special divisor Pi + ■ ■ ■ + P g +i — oo + — oo_ is not 
linearly equivalent to a positive divisor. 

Then, for \s\, \t\ sufficiently small, the linear vector space of Baker- Akhiezer eigenf unctions 
associated with (E,P, oo±) is 2- dimensional and has a unique basis ip 1 and if) 2 , such that 



the vector-valued function i/j = 




has the following normalized expansions at oo_ and 



oo +; respectively: 

if) = e -iO+2A 2 t) 

Remarks: 

1. The Genericity Condition can be rephrased as, there exists no non-constant meromor- 
phic function with pole divisor P\ + ... + P g +i which vanishes simultaneously at oo + and 

at oo_. 

2. The proof of uniqueness follows from the fact that the ratio tp 1 /^ 2 of two Baker- 
Akhiezer functions is a meromorphic function whose poles lie in the zero divisor of ip 2 . The 
condition of non-speciality ensures that the dimension of the space of such functions is 2. 
Since l(T> — oo + ) = l(T> — oo_) = 1, we can normalize ip l and tp 2 by requiring that they 
vanish at oo + and oo_, respectively. Because l(V — oo + — oo_) = 0, which follows from the 
Genericity Condition, they must be independent. For more details, see £ 15.21 

The existence of the Baker-Akhiezer eigenfunction is shown by explicitly constructing 
its components in terms of the Riemann theta function of S. We summarize the main 
steps of this algebro-geometric construction. (See [TJ H2] for a comprehensive discussion and 
application of this method to several other soliton equations.) Let ai, ... ,a g , b±, ... , b„, be 
a homology basis for S, such that a« • Oj = 0, 6j ■ bj = 0, Oj • bj = 5y (see Figured)). 
Let u>i, ... u g be holomorphic differentials on S, normalized as follows 

Uj = 2ni5 jk , j, k = 1, g. 

We introduce the period matrix B of entries 

B jk = <b Uj, j,k = 1, ... ,g, (9) 

Jh 

and construct the associated Riemann theta function 

0( z ) = ex p (K n ' Bn ) + ( n ' z » . ze £9 - ( 10 ) 

n<EZ9 



+ 0(A 



-1\ 



+ 0(A 



-1> 



n 
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Figure 1. Homology cycles on £, adapted from |T]. The dashed vertical 
lines are branch cuts, the solid curves are portions of homology cycles on the 
upper sheet, the dashed curves are portions on the lower sheet, and the dotted 
horizontal line is the real axis in the A-plane. Also shown are cycles a g+ \ and 
£, used in §|Uand|3J 



The series ()10|) is absolutely convergent; this follows from the fact that Re(£>) is a negative 
definite matrix. Moreover, if e^'s are the standard basis vectors of C 9 and = Bek, for 
k — 1, . . .g, then 

0(z + 27rie fe ) = 0(z), 0(z + f fc ) = e-^ Bkk - Zk 6(z); 

in other words, 9 is a quasi-periodic function. (The e^'s are called periods and the f^'s the 
quasiperiods of 6.) 

We also introduce the Abel map, which takes a divisor Y2k=i Qk to a point of the complex 
torus Jac(S) = C 9 /A, where A is the 2g-dimensional lattice spanned by the columns of the 
matrix (I \B). The Abel map of £ is defined by 

_£_ 9 rQk 

^■Z2^ k ^z2 ^modA, (11) 
fc=i k=i J p ° 

where Pq G S is a given base point. 



2.1. Eigenfunction formulas. The components of the vector- valued Baker- Akhiezer eigen- 
function i/> are constructed in terms of ratios of Riemann theta functions: 

M s,f) .e Xp ( is(Sil( P) - f ) + TO + f )) ^fl^±^M 

0(A(oo_) - - K) g + {P) 



0(A(oo_) + iVs + iWt - A(V + ) - K) #+(oo_)' 

(12) 
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and 

M s,t) = e*p {W) + f ) + IMP) - f )) e{A(P l^) + -*V:*V~ K) 

9(A(oo + ) - A(V-) - K) g_(P) 



9(A(oo + ) + iVs + iWt-A(V_)-K) g-(oo+) 

(13) 

The ingredients of this formulas are described below: 

(1) The essential singularities of if) at oo± are introduced by means of the unique Abelian 
differentials dQi(P) and df^C-P), defined by their asymptotic expansions at oo± 

dfii ~ ±dA, dfi 2 = ±4AdA, as P -> oo±, 

and normalized so that their integrals vanish along the a-cycles. 

(2) We will choose Pq = oo_ as basepoint for the Abel map A, and choose the rightmost 
branch point \2 g +2 m the lower half-plane as basepoint for the Abelian integrals 



Qj(P) = / d£lj. This latter choice gives the property 

J A2 s +2 



Qj o i = —Qj, 

where t : (A, /i) i— > (A, — /i) is the sheet interchange automorphism. In this way, the 
asymptotic expansions for flj and 1^2 at oo± are 

fii = ±(A-f + o(l)) and fi 2 = ±(2A 2 + f + o(l)), 

for certain constants and AT. 

(3) Although they have different basepoints, we make the convention that the paths 
of integration for A(P) and Qj{P) are yoked, so that a homology cycle may be 
added to one path only if it is added to the other; in other words, the difference 
between the paths for Qj{P) and A{P) is a fixed path from oo_ to \2 g +2- The 
most straightforward way of choosing this is to cut the Riemann surface along the 
homology cycles aj, bj, resulting in a cut surface S with boundary aj + bj — aj—bj, 

j 

and fixing a path II between oo_ and A 29 +2 in the interior of S (see Figure |2J). 

(4) Then frequency vectors V and W are chosen to make ip 1 ,^ 2 well-defined functions 
on S. In fact, if the paths of integration are modified by adding an integer linear 
combination XwUi n k<ik + °mkbk of homology cycles, then xj> changes by the factor 



exp 



E 



rrik ( is <p dQ\ + it (b dfi 2 ) — ™>k (isV^ + itWk 

b k Jb k 



,fc=i 

which equals 1 if we define the components of V and W as 



V k = <b dflx, W k = <b dn 2 . 

Jb k Jb k 

(5) T>± are the unique positive divisors linearly equivalent to V — oo±, and K is the 
vector of Riemann constants ^2] which has the property that the zero divisor of 
6{A{P) - A{V ± ) - K) is V±. 
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(6) In order to make ip 1 and ip 2 have pole divisor T>, the components of ip are multiplied 
by meromorphic functions g± whose zero divisors are V± + oo T and whose poles lie 
in the original divisor T>. For the sake of definiteness, we will normalize 

9+(°°-) = 9-(°°+) = !■ 




Z 





2g+2 




Figure 2. Cutting a two-hole torus E along a homology basis, to yield the 
cut surface S 

For later use (e.g., deriving formulas for the components of the VFE solution curve 7 given 
by ©), we will explain how to derive a more compact and convenient expression for the 
eigenfunction. Let 

roc + 

D = A{V + ) + 7C, r = / u, 



both computed using paths in So. Then A{T>J) + 7C = D + r, so 

? = 9+{P) ex P H n l{P) - f ) + TO + f )) ' W> + lV * + ™ - D > ^ 

^ 2 = g_(P) exp (is^P) + f ) + it(Q 2 (P) - f )) 



8(iVs + iWt - D) 6(A(P) - D) ' 
e{A{P) + iVs + iWt - D - r) 0(D) 



8(iVs + iWt - D) 9(A(P) — D — r) 
Next, introduce another meromorphic differential cir2 3 , uniquely defined by 



an 3 ~ ±^ 

A 



as P — > oo±, 



and the requirement that it has zero a-periods. Let ^(-P) be the corresponding Abelian 
integral on E, computed with the same basepoint and convention for paths as fl\,fl 2 . Then, 
by imposing reality conditions that guarantee that the VFE solution is a curve in R 3 at all 
times (see ^5. 21 for details), we obtain the following simplified expression for the components 
of ij): 



^ = exp (is(lMP) - f ) + it{n 2 (p) + f )) 
ij 2 = -iexp (is(0!(P) + f ) + it(Q 2 (P) - f ) + Q 3 (P)) 



0{A{P) + iVs + iWt - D) 



6(iVs + iWt - D) 

9(A(P) + iVs + iWt — D — r 



(14) 



6(i~Vs + iWt - D) 



(15) 
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An explicit computation involving the asymptotic expansions of the eigenfunction if) at oo± 
(see, e.g., |J) shows that ij) is a common solution of the pair of NLS linear systems (j3J) with 
potential q given by 



q{x,t) = Aexp(-iEs + iNt) 
where the formula for the constant A is given in 



9(iVs + iWt - D + r) 
6(iVs + iWt - D) 



(16) 



2.2. Curve Formulas. When 7r(P) is real, a fundamental solution matrix of (jlj) is given 
by 



^(P) -jj 2 (P ) 
ip 2 (P) ^(P) 



(17) 



The Sym-Pohlmeyer reconstruction formula (jSJ) yields the following expression for the com- 
ponents of the skew-hermit ian matrix T representing the position vector of the VFE solution: 



r n = i(n{p)s + n' 2 (p)t) 
i 



r 



2det# 



l^l'viog ( 9{A{P) + ^) + i^Tviog ( d{A{P) + y " r) 

W 1 V ° g V(9(^(P) - <p)J + IV> 1 V ° g ^(A(P) - V - r) 



dA{P) 
dX 



21 



det^V V M-OH + V) ; dA 31 ; 



(18) 

where y> = iVs + iWt - D, and Q'^P) = dfy(P) /dX. 

It is easy to verify that, given a fundamental solution matrix \l/ of the linear systems (@J), 
the quantity det \I/ = \ip x \ 2 + \ip 2 \ 2 is independent of both s and £. Differentiating it with 
respect to s, one gets 



(6(A(P) + y) \ /gGjcg + y - r) 

1 V g V 0(.4(P) - V ; + IV 1 V ° g v 0(A(P) - <p - r) 



2 det W log0(y>) 



■V = 0. 



Differentiation with respect to t yields a similar identity with V replaced by W. Since the 
two frequency vectors span C 2 , one concludes that the quantity between the square brackets 
must vanish and obtains the following 

Proposition 2.3. The components (71, 72, 73) of the position vector*/ of an N-phase solution 
of the Vortex Filament Equation are given by the following expressions 



21 



7i + 172 



ip 1 ^ 2 f (0{A{P) +iVs + iWt-D - 
det¥ V ° g V 0{A(P) + i\s + iWt-B) 



r n = 73 = i[n[(p) s + n' 2 {p)t] + viog#(iv s + iwt - d) 

where vr(P) = A G R is the reconstruction point. 



dA{P) 
dX 
dA{P) 

dX ' 



+ n' 3 (P) 



2.3. Closure Conditions. The potential q associated with a closed curve of length L by 
the Hasimoto map is L-periodic, but only up to multiplication by a unit modulus constant: 



q{s + L)= e^q(s) 



where 



t ds. 
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However, a quasiperiodic potential like this is related to a genuinely periodic NLS solution 
by the NLS symmetry mentioned in the introduction: 

q{s,t) =e l(as - a2t) q{s-2at,t). 

(Obviously, one chooses a = —<fi/L to make q a periodic potential.) This symmetry lifts to 
the NLS linear system: 

t{)(s, t; q, A) = exp(|(as — a 2 t)a 3 )ip(s — 2at, t; q, A), where A = A — |a. 

(This transformation, for fixed time t = 0, appears in Grinevich and Schmidt [16J.) It follows 
that the continuous spectrum of q is the same as that of q, but translated by —a/2, and the 
Sym formula © yields the same curves at corresponding values of A: 

rU=A = r|A=A -|a- 

Thus, we can without loss of generality assume that a closed curve of length L is generated 
from a potential q(s, t) that is L-periodic in the variable s. Note that, for finite-gap solutions, 
the value of the frequency vector V is unchanged under horizontal translation Xj i— > Xj — |a of 
the branch points, but the phase constant E changes by E i— > E — a. In fact, our assumption 
that q is L-periodic implies that E has been translated to zero and the entries of V are 
rational multiples of 2tt/L. Therefore, if V is the greatest common divisor of the entries of 
V, then q has period T = 2tt/V, and L is an integer multiple of T. 

Since the natural curvatures (which are one half the real and imaginary parts of q(s,t)) 
are periodic, the corresponding natural frame vectors must also be periodic. Since those 
are obtained by conjugating a fixed basis of su(2) by the fundamental matrix it follows 
that \I/ must be L-periodic or antiperiodic. For finite-gap solutions satisfying the periodicity 
condition, this amounts to assuming that Q\(P) is an integer multiple of tt/L. 

Examining the formulas (|18|) shows that our assumptions so far guarantee that r 2 i is 
periodic, and that r n is periodic if and only if dfii(P) vanishes. We summarize these 
closure conditions as: 

Proposition 2.4 (Closure). Suppose q(s,t) is a finite-gap potential of period 2ir/V in s. 
Then the curve reconstructed from q using the Sym formula (JUJ) with A = X(P) is smoothly 
closed of length L = 2nrr /V if and only if (a) dQi(P) = and (b) exp(iLf2i(P)) = ±1. 

The points which fulfil the closure condition (b) are discrete points of the Floquet spectrum 
associated with the potential q, which can be described in terms of Floquet discriminant. 
The Floquet discriminant of q is defined as the trace of the transfer matrix of ^ over the 
interval [0, L\: 

A(q; A) = tr (^{s, t; A)~ 1 \I>(s + L, t; A)) . 

(Because the transfer matrix only changes by conjugation when we shift in s or t, A is 
independent of those variables.) Then, the Floquet spectrum is defined as the region 

a(q) = {X E C | A{q; A) G E, -2 < A < 2}. 

Points of the continuous spectrum of q are those for which the eigenvalues of the transfer 
matrix have unit modulus, and therefore A(q; X) is real and between 2 and —2; in particular, 
the real line is part of the continuous spectrum. Points of the L-periodic discrete spectrum of 
q are those for which the eigenvalues of the transfer matrix are ±1, equivalently A(q; A) = ±2. 
Points of the discrete spectrum which are embedded in a continuous band of spectrum have 
to be critical points for the Floquet discriminant (i.e., dA/dA must vanish at such points); 
thus, the conditions in Proposition 12.41 mav be rephrased as requiring Aq to be a real zero of 
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the quasimomentum differential that is also a double point. These conditions coincide with 
those derived by Grinevich and Schmidt |16j using different reasoning. 

For finite-gap potentials, we may use the simplified eigenfunction formulas p4|1 . (jl5j1 to 
write the Floquet discriminant as 

A(q; A) = 2cos(PT2 1 (P)), A = tt(P). 

If we have an explicit formula for the discriminant, we can use it to locate the points of 
the discrete spectrum; see §3.51 for an example. However, it is very difficult, in general, to 
determine how to choose the branch points so that a point of the discrete spectrum coincides 
with a zero of the quasimomentum differential. Even in genus one, this corresponds to an 
implicit equation involving elliptic integrals. 

2.4. Periodic Self-Intersection. Whether or not the filament closes smoothly, it may 
intersect itself after one period T of the potential. (For example, this happens with elastic 
rod centerlines [20j.) In this section, we will calculate the general condition under which 
this happens. In the expression for Tn in (|18|). the only term which is not automatically 
T-periodic is the linear term iQ[(P)s. Therefore, Tn is T-periodic if and only if f^(P) = 0, 
i.e., the first condition of Prop. 12.41 above is satisfied. So, we will assume this condition. 
Next, r 2 i is T-periodic if and only if 

«V = -iexp(2 l( n 1 (P). s + n 2{P )t) + n 3{ p)) 0{A(F) + v) !^ P) + v - t) 

6(<p) 2 

is periodic. But that amounts to assuming that A(P) is a point for the T-periodic discrete 
spectrum, giving a smoothly closed filament of length T. However, self-intersection may also 
occur if the last factor in r 21 vanishes for some s-value: 

f9(A(P) + <p-r)\ dA{P) 



( 6(A(P) + y-r) \ dA(l 
g V e(A(P) + y>) J dX 



Q' 3 (P) = 0. (19) 



(Note that this factor is also T-periodic in s.) Then r 21 is not T-periodic, but is equal to 
zero at a succession of s-values, and the T-periodicity of Tu guarantees that the filament 
returns to the same point in R 3 . 

For finite-gap solutions, the first self- intersection condition Q[(P) = is polynomial in 
A. In general, the second condition (JT§j) depends on A, so and also t, so that this kind of 
self-intersection may only be present at one time. However, in genus one, a shift in time 
may be offset by a shift in sq, so that these self-intersections persist in time. Furthermore, 
in genus one the derivatives of the logarithms may be expressed in terms of Jacobi zeta 
functions, and we may choose s and t to give cp = 0. This yields a polynomial condition in 
A. We conjecture that an analogous polynomial condition occurs in higher genus. 

3. Genus One Solutions 

In this section, we will work out the representation of solutions of the filament flow in 
detail for the case when the genus is one, when E is an elliptic curve given by the equation 

^ = n ( a - ^ = ( a - Ai )( a - a2 )( a - Ai )( a - A2 )- ( 2 °) 

3=1 

The Floquet spectrum of a generic 2-phase solution is shown in Figure EJa). In particular, 
the intersection points a.\ and a 2 of the complex bands of spectrum with the real axis are 
the real zeros of the quasimomentum differential dQi given in equations (J3H)) . (J3~%j) below. As 
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discussed earlier, closed curves are obtained by selecting the branch points Xj so that one of 
the aij is also a double point, and using A = <x, in the reconstruction formula (JHJ). 

3.1. Elastic Rods. Filaments generated by genus one potentials are interesting objects 
from geometric and topological points of view, because they are the centerlines for Kirchhoff 
elastic rods. These are thin rods that are critical for an energy that consists of one term for 
the isotropic bending energy (i.e., a multiple of J K 2 ds, where k is the Frenet curvature of 
the centerline) and a term for the twisting of the rod about its centerline. Solutions of this 
physical variational problem are exactly the solutions of a geometric variational problem, 
where the Lagrangian is of the form 



(The values of the parameters h, I2 depending both on the relative weight of the bending and 
twisting energies and on the boundary conditions in the physical variational problem; see 
[2*U] and [23 for details.) Regarding /1, 12 as Lagrange multipliers, we can think of Kirchhoff 
elastic rod centerlines as extremizing the bending energy subject to constant total torsion 
and constant length constraints. These curves include classical Eulerian elastic curves (for 
li = 0) and free elastica (for li = l 2 = 0). 

The identification between travelling wave solutions of the NLS and elastic rod centerlines 
was first discovered by Kida j22], who sought solutions of the filament flow that moved by 
rigid motion (i.e., rotation and/or translation) plus possibly a translation along the curve. 
(As we will see below, under the Hasimoto map these filaments correspond exactly to finite- 
gap NLS potentials of genus one.) Kida obtained parametrizations for such curves, in terms 
of elliptic integrals, showed that there exist smooth closed curves of this type, and identified 
the planar curves of this type with Eulerian elastica. 

Later, Langer and Singer [27] began with the variational problem for T above, and showed 
that solutions were precisely those curves which move by rigid motion (plus translation along 
the curve) under the filament flow. Langer and Singer expressed the curvature and torsion 
of these curves in terms of periodic elliptic functions; consequently, these curves consist of 
repeated congruent segments which smoothly close up if certain elliptic integrals satisfy a 
rationality condition. Figure El shows some examples. 



Like these specimens, all closed elastic rod centerlines have discrete rotational symmetry 
about an axis. It is not hard to see that if such a curve is knotted, it must be a torus knot. 




(21) 




Figure 3. Some examples of knotted elastic rods. 
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Ivey and Singer [20] later showed that every torus knot type is realized by a one-parameter 
family of smooth elastic rod centerlines. In fact, each such family connects smoothly to 
another family of centerlines with the same rotational symmetry, but with a complementary 
torus knot type (like (2,5) and (—3,5) shown above). We call this a homotopy of closed 
elastic rod centerlines. 

3.2. Genus One Potentials. We now turn to the detailed calculation of genus one NLS 
potentials. The data necessary to construct q will be expressed in terms of elliptic integrals 
and algebraic functions of the branch points. For these calculations, we rely heavily on 
the treatise by Byrd and Friedman jlj; citations of the form "formula xxx.yy" refer to this 
volume. 

We begin by choosing a standard homology basis a, b on the elliptic curve £ (see Figure 
0J) and computing the periods of holomorphic differentials. 




Figure 4. Branch points for elliptic curve E with oriented homology cycles 
and branch cuts. The upper sheet tends to oo + to the right and left. 



Lemma 3.1. 



la ^ |Ai — A2I Jb ^ I Ai — A2 

where we introduce the abbreviations 

K := K{p), K' := K{j,') 
for complete elliptic integrals of the first kind with moduli 



(22) 



lAi-A 



p ' := nT — T~\' P-=V 1 -(p') 2 - 
|Ai — A2I 

Proof. We will give a somewhat detailed explanation of this calculation, since it is a model 
for others below. 

To transform the integrals in (|22|) to elliptic integrals along the real axis, we use the linear 
fractional substitution 

A(*) = A 2 + ^_^, h:=*-^, (23) 
1 nx — A 2 

under which Ai, Ai, A2, A2 correspond to x-values 1, {p , )~ 2 , 0, 00 respectively. Then J2QJ) gives 

2 

f = -h 2 \X, - A 2 | 2 (yf^) (x{x - 1)(1 - (p'fx)) ■ (24) 
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Without changing the value of the a-integral in ([22)1 . the a-cycle may be contracted in the 
A-plane to a straight-line path that travels on the left-hand edge of the branch cut, from Ai 
to Ai, and then back along the right-hand edge of the branch cut. Since /i differs only by 
a minus sign on the two sides of the branch cut, the a-integral is twice the integral along 
the left-hand edge, where Re(/i) > 0. As A moves along this path, x moves from (p')~ 2 to 1. 
Thus Im(dA) > and Re(dy) < as we move along this path, so the equation 

1 - hy 

shows that Re(/i/(l — hy)) < there. Hence, along this path we take 

-h 



V = |(Ai - A 2 )(A 2 - A 2 )| [y^) " " (pO 2 *), 

giving 

1 fdX _ -2i f 1 dy 

2Ja V ~ [Ai — A 2 | J {pT 2 y/x(x-l)(l-(j/) 2 x)' 

Reversing the limits in the integral on the right-hand side yields the complete elliptic integral 
K(p) (see formula 236.00). 

The computation of the o-integral is similar. □ 

Consequently, the normalized holomorphic differential (with a-period equal to 27ri) is 

„ . /|Ai-A 2 |\ dA 

and the 'Riemann matrix' for the ^-function used in the construction of q is 



B = ju = -2nK'/K, (25) 
The frequencies inside the ^-function are 

V = 4vri ( {dX/fj) =7i\X 1 -X 2 \/K 



and W = cV, where 

2g+2 

c= X i = Ai+Ai + A 2 + A 2 . (26) 

Using the formulas (j64|) . ()65jl . the frequencies in the exponential part of q are given by 

E = c-2L, N = 2cL-Ad-c 2 , 
where c is as in ()26|) . d is given by (jHSJ), and L is the following ratio of a-periods: 

JaV I J a V 

To calculate the numerator of L, we use the same path and change of variables as in the 
proof of Lemma f3. II Thus, along the left edge of the branch cut, we have 

^dA= =L (x 2 + X2 - J ^ dy 



V ' |Ai — A a | V 1-hxJ y/ x (x - 1)(1 - (p') 2 x) ' 



FINITE-GAP SOLUTIONS OF THE VORTEX FILAMENT EQUATION, I 15 

Applying formula 236.02, we get 

L = A 2 + (A 1 -A 2 )H^), (27) 

where U((3 2 ,p) denotes a complete elliptic integral of the third kind, with parameter 

2 ■= p2 = Al ~^ 
1-h Ai-A 2 ' 

Note that we may also calculate L using formula 235.02, obtaining 

L = A 2 + (A 1 -A 2 )^%^. (28) 

This is equivalent to (|2*7|) by application of the addition formula 117.03 for elliptic integrals 
of the third kind. Furthermore, comparing (|27jl and ()28jl shows that L is real, as it must be 
since complex conjugation 

r : (A,//) i (A, /I) 

takes cycle a to —a. 

To calculate the shift r in (|16p. we use the Abelian differential 

dfi 3 = ^— ^dA 
H 

which has a-period equal to zero. Using formula 233.02 to compute £(A///)dA gives 

r = - hn 3 = 4 [(A 2 - L)A" + (A 2 - A 2 )n(/i,p')] • (29) 
A |Ai — A 2 | 

One can verify that lm(r) = 7r by using the addition formula 117.02 to simplify U(h,p) + 
H(h,p), given that h = (p') 2 /h. 

It will later be convenient to have an alternate formula for r. Using addition formula 
117.05, we can write 



;i - h)K (U(h lP ') - K') = hK'Tl(p 2 ,p) - \\\y\F{^'), = arcsin (Vh/p') 



where F denotes an incomplete elliptic integral of the first kind, e.g., F(ir/2,p') = K' . (Note 
that \h\ = (p') 2 < 1, arg(/i) G [0, 7r), arg ((1 - h)/(h - 1)) = 2 arg(A 2 - Ax) G (0, vr), and 
we take the square root to have positive real part.) We use this formula to substitute for 
U(h,p') in (}2T?j) . and use (|2*5jl to substitute for L, to get 

r = -F^p') (30) 

Finally, the modulus of scalar A may be calculated by numerical integration as described 
in Appendix 15.31 However, we also have the following formula, 



,,, |Ai-A 2 | 2pp'K 

|A| = ^rv^' (31) 

which may be derived by substituting the formula (jlfij) for the potential into the NLS and 
evaluating at s = t = 0. 
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3.3. Geometric Interpretation. The above information is enough to enable us to discuss 
the geometric invariants (i.e., curvature and torsion) for the filaments corresponding to a 
genus one potential q(s,t). 

Because a translation in time modifies q only by a unit modulus factor which is independent 
of s, and may be absorbed inside the theta function by translation in s, the geometric 
invariants at time t differ from those at time zero only by a constant-speed translation in s. 
Thus, the filament at time t is obtained from the filament at time zero by Euclidean motions 
(accompanied by a translation along the curve). It follows from the semigroup property of 
the evolution equation that the Euclidean motions must belong to a one-parameter subgroup. 
Together with the discussion at the beginning of this section, this gives the following 

Proposition 3.2. All finite-gap solutions of genus one for the vortex filament equation move 
by rigid motions, and are congruent to elastic rod centerlines. 

It is instructive to express the curvature of genus one solutions in terms of elliptic functions. 
By the above argument, we may assume that t — without loss of generality. Also, since 
the reality condition (see Appendix 15 .2|) implies that D is purely imaginary, we may assume 
that D = 0. 

We first observe that for B given by (J2HJ), 

e(2iz) = e 3 (z) = e (z + 7r/2), 

where df. denotes the Jacobi theta function with modulus p and period 2ir (for k — 1,2) or 
7T (for k = 0,3). Next, let p = |Re(r). Then 

For convenience, introduced the rescaled variable z = (Vx+n)/2. Then, using the translation 
properties and addition formulas 1051. ff for Jacobi theta functions, we have 

2 ro^-nM2 ln A Mz + i P )6 3 (z - ip) . A |2 0f(ip) / 0f(i P )0 2 (*) 



kT = |2g(s,0)r = \2A\ Z oy ^ = \2A 



el{z) 1 1 £ 2 (0) V eiii P )ei{z 

This formula coincides with the curvature formula for elastic rod centerlines derived in 
for, if we set 

w := dc(2Kip/7r) = dn(2Kp/n,p'), (33) 
and substitute in for \A\ from (|3*Tj). then we have 



k 2 = ( 1 — sn 2 u) , where u = z = — (s + n/V), 

1 W Z I 7T 2w 



and we use /to = 2w|Ai — A 2 | to denote the maximum curvature. 

We will now make some specific observations about the relation between the geometry of 
the filament and the spectrum in genus one. For these, we will need the following formula: 



Lemma 3.3. 



w2 = IMA, -3)|' 
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Proof. Using (|30|). we compute that 

w = dn{2Kp/7r,p') = ics(Kr/ir,p') = icn(2F,p')/ sn(2F,p') 
where F = F(<f),p'), as in (jSOjl . Then, using double angle formulas 124.01, 

sn 2 (F,p') dn 2 (F,p') - cn 2 (F,p') Re(h - 1) 
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2isn(F,p / )cn(F,p')dn(F,p / ) \h-l\ ' 

where the ambiguity of sign arises from solving sn(F,p') = Vh/p' for the values of cn(F,p') 
and da(F,p'). Then substituting (|2~3J) yields the result. □ 

In |2Zj, it is shown that the curvature and torsion of an elastic rod centerline satisfy 

2r-h=j/K 2 , where j 2 = - w 2 )(w 2 - p 2 ), (35) 

and l\ is the ratio of the coefficient of the term j r ds to the coefficient of the elastic energy 
J jK 2 ds in the Lagrangian ([21)1 . It is clear that the torsion of the centerline is constant if and 
only if w 2 = p 2 or w 2 = 1. From (J34j) we see that these correspond to the cases where Ai — A2 
is purely real or purely imaginary, respectively. In these cases, we may translate the branch 
points so that they are symmetric under the reflection A 1— > —A about the origin. (Recall 
that the curve produced by the Sym-Pohlmeyer formula is unchanged provided we apply the 
same translation to the reconstruction point Ao.) We summarize these observations: 

Proposition 3.4. A filament associated to a finite-gap solution of genus one has constant 
torsion if and only the branch points are symmetric about an axis perpendicular to the real 
axis. 

In §3] of this paper, we will discuss what symmetry of the branch points means in higher 
genus. 

3.4. Quasimomentum Differential. Before making further observations, we need to cal- 
culate the meromorphic differential dOi for genus one. (Recall that the zeros of this 
differential are used to reconstruct closed curves.) This will later enable us to calculate the 
Floquet discriminant, and (in Prop. 13.51 below) to connect the coalescence of the a^s with 
the filament being an elastic curve. 
The form of this differential is 

dOx = A " ~ f A ~ c i dA> (36) 

where c is as in (J26|) and c\ is chosen so that the a-period of dfii is zero. The value of 
§ (A//x)dA was obtained in the above calculation of L. In a similar way, we apply formulas 
236.17 and 336.02 to get 



-dA 41 



a V Ai - A 



2 



A 2 K + 2A 2 (A 1 -A 2 )n(/? 2 ,p) 



+ 2(1 l A ^f p2) (f^P) + (P 2 ~ P) K + (20V + 2/5 2 - /3 4 - Sp 2 M(3 2 , P )) 



(37) 



where E(p) is the complete elliptic integral of the second kind. We combine the previous 
two results to compute 

Cl ./^ dA //«. 

Ja V I Ja V 
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Miraculously, the coefficient of IT(/3 2 ,p) in C\ is zero, giving 

" dA 



A — | A — \ ( |Ai — A2 1 2 ^ — |Ai| 2 — [ A2 



(38) 



The roots of the quadratic polynomial in A enclosed in square brackets are the zeros a\, 0*2 
of the quasimomentum differential. 

Proposition 3.5. Suppose a genus one filament is a quasiperiodic space curve, i.e., it has 
periodic curvature and torsion, and is bounded in space. Then the filament is an elastic curve 
if and only if ai = 02- 

Proof. In j^Zj (see equation (26) in that paper), it is shown that elastic rod centerlines are 
bounded in space if and only if 

mv) . 2 a.t^zf 

__ 1 = W -p + — g_ 

K K 2 Q 

where l\ is the coefficient in the geometric Lagrangian (|21j) for elastic rod centerlines. Thus, 
quasiperiodic elastic rod centerlines are elastic curves if and only if 

__ 1 = w - p (39) 

On the other hand, the discriminant of the quadratic in (J??%|) is 

j + 4ci = J(Ai + Ax + A 2 + A 2 ) 2 + 2|Ai - A 2 | 2 ^ - 2|A a | 2 - 2|A 2 | 2 . 

By using (jMj) . one sees that this vanishes if and only if (|59~j) holds. □ 

Remark 3.6. The assumption of quasiperiodicity seems to be necessary here. For, without 
changing the spectrum of the curve, we can shift the reconstruction point Ao in the Sym- 
Pohlmeyer formula; this will add a constant to the torsion, causing the curve to no longer 
be quasiperiodic, and also causing the curve to cease to be an elastic curve. 

3.5. Floquet Discriminant. Following our earlier discussion the trace of the transfer ma- 
trix over one period of q will be given by 

A (A) = 2 cos (yQ^P) 

where P is a point on the hyperelliptic curve lying over A and Oj(P) is the integral of dfli 
from basepoint A2 to P. This integral is path-dependent, but since dfl\ has zero a-period 
and 6-period equal to V, the right-hand side is unaffected by the choice of path. It is also 
unaffected by the choice of branch for P, since i*dfli = — dfli, where t is the sheet interchange 
involution. 

Because the integral formulas used in §3.21 are better suited to using basepoint Ai, we 
will do that instead. Because fj 1 dfli = jj^ dQi = and dfii = ±V/2, this change of 
basepoint changes the value of A(A) by a minus sign. 

Let r denote the path chosen from Ai to P. Initially, we will assume that P and T lie on 
the upper sheet, to the left of the branch cut from Ai to Ai. (On this sheet, Re(/i) > above 
the real axis.) However, the formula for A(A) which we will derive extends analytically for 
all A. 
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Let X denote the x-v&hie corresponding to A under the change of variables. Then by a 
calculation similar to the proof of Lemma 13. II (in particular, using 236.00 in [3]), 



dA 
r V 



|Ai — As 



dy 



2i 



(p'Y 



^x{x-l){l-{p>) 2 x) \\i - A 2 



where 



ip = arcsin I y 1 — X(p') 2 /p 



Let u = F(ip,p). Note that this is an analytic function of if) provided |Re(sin-^)| < 1/p. 
In fact, this is true if we restrict X to lie strictly inside the circle in the x-plane centered on 
the real axis and passing through 1 and (p')~ 2 - Then we may take 



sn« 



y/l-X(j/)y P , 



cnu 



p 



dnu = p'y/X. 



(40) 



By a calculation similar to that yielding L, we get 



-dA 



2i 



At — Ao 



r H 



I A 



where U(ip , (3 2 , p) denotes an incomplete elliptic integral of the third kind. However, in cal- 
culating the T-integral of the quasimomentum differential dfli, the coefficients of H(ip, (3 2 ,p) 
cancel out, so we will omit these terms from now on. 
Applying formulas 236.17 and 336.02 again, we get 



V 2 

-dA 

r fJ> 

2i 



Ai — A2 



\r,U + 



(Ai-A 2 ) ; 



2(1 - (32)^2 _ p2) 

Finally, combining these terms gives 

. A 2 

dOi 



f3 2 E(ij,p) + (p 2 -p 2 )u-p 



cn u dn ksiik 
1 — f3 2 sn 2 u 



■ (41) 



)\ - Ci 



dA 



i|Ax-A 2 | [Z{u)-0 



, cn u dn u sn u 
1 — ft 2 sn 2 u 

where Z{u) denotes the Jacobi zeta function with modulus p. Consequently, the Floquet 
discriminant is 



A(A) 



-2 cos 



2\K Z(u) - 0- 



, cn u dn u sn u 
1 — ft 2 sn 2 u 



1 1 



F ^arcsin ^y 7 ! — X{p') 2 jp 



(42) 

(The minus sign in front arises from our change of basepoint from A 2 to Ai.) 

In an earlier paper 0, we made a similar calculation of the discriminant for genus one 
solutions of the sine-Gordon equation. Like that formula, the argument of the cosine in A(A) 
here consists of a Jacobi zeta function minus an algebraic term. For, by substituting in from 
(l4*U|) and comparing with (j23), we obtain 



0- 



. cn u dn u sn u 



;i-/3 2 sn 2 w) IA1-A2I VA-A 
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where sign may be checked by verifying that the left-hand side has positive real part when, 
say, X = 1/p'. (Recall that we are assuming that Re(/x) > 0.) For the purposes of in- 
vestigating the discriminant numerically, this alternate expression is not useful because the 
branch cut involved in calculating ip does not necessarily coincide with the branch cut for 
/i. However, the simultaneous branch cuts in the two terms in (|42J1 cancel out, modulo 2tt, 
showing that A is a smooth function of A. 

The spectrum may be plotted using level curves of Re(A(A) and lm(A(A)); for, the spec- 
trum is contained in the level curve lm(A(A)) = 0, the simple points (which lie at the end of 
the continuous spectrum) lie at the intersection of this level curve with a level curve where 
Re(A(A) = ±2, and points of order k lie at the intersection of k such level curves with the 
continuous spectrum. 

About Figured This figure shows the Floquet spectra of some elastic rods (including some 
with special geometric features) belonging to a typical homotopy. The diagrams, which are 
based on Mathematica plots using equation (J4*2~j) . show the continuous spectrum, including 
the real axis, the simple points (which appear as dots), the multiple points (which are usually 
double points, and appear as x's), and the Sym-Pohlmeyer reconstruction point (which is 
circled). 1 In each case, the spectrum is translated in the real direction so that the sum of 
the branch points is zero. The stages of the homotopy shown are: 

(a) an unknotted rod, with branch points — 0.454 ±0.095i and 0.454 ±0.324i in the upper 
half plane; 

(b) a constant torsion rod, with symmetric spectrum and branch points ±0.388 ± 0.316i 
(note the coalescence, at the base of the left-hand spine, of two double points, one of 
which previously entered from the left); 

(c) another unknotted rod, with branch points -0.342 ± 0.409846i and 0.342 ± 0.320i, 
and double points migrating out along the left-hand spine; 

(d) a third unknotted rod, with branch points —0.206 ± 0.545i and 0.205 ± 0.366i, and a 
double point which has entered from the left; 

(e) an Eulerian elastic curve, with branch points — 0.201 ±0.546i and 0.201 ± 0.370i, and 
a multiple point of order six at the origin, where the zeros of the quasimomentum 
differential have coalesced; 

(f) a knotted rod, with branch points —0.251 ± 0.557i and 0.251 ± 0.308i, and the recon- 
struction point now on the left-hand spine. 

As the homotopy continues, the right-hand spine shrinks toward the real axis. After the 
change in knot type, the spectrum does not look qualitatively different from frame (f). 



When the curve is smoothly closed, this point is necessarily a multiple point of order at least four |10| . 




Figure 5. Drawings of the spectra of selected closed elastic rods belonging to 
a homotopy between a singly- covered circle and a double-covered circle, going 
from unknotted to a trefoil knot type along the way. 
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Figure 6. Drawings of the the rod centerlines corresponding to the spectra 
in Figure El 



finite-gap solutions of the vortex filament equation, i 
4. Symmetric Solutions 
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As before, let £ be a smooth hyperelliptic curve of genus g, defined by 

2g+2 

J'=l 

where the branch points are symmetric about the real axis, i.e., X g +k = for 1 < k < g+l. 
In this section, we assume that the branch points are also symmetric about the origin. Then, 
in addition to the sheet interchange automorphism t, £ has an automorphism 

a:(A,/z)~(-A,(-l)* + V), (43) 

the sign being chosen so that a fixes the points oo + and oo_ on £. We shall refer to such 
hyperelliptic curves, and the NLS potentials that arise from them, as symmetric. 

4.1. Consequences of Symmetry. In order to see how the standard homology basis be- 
haves under automorphism a, it is convenient to define an extra cycle a g+ i, as shown in 
Figure [TJ Then, in terms of the a,j for j < g and the clockwise cycle I about oo + , 

9 

a g+ x ~ -I- J2 a J- ( 44 ) 



3=1 



(Here and below, ~ indicates equivalence of basepoint-free homology cycles within £\{oo + , oo_ 

Proposition 4.1. For 1 < j < g, cr(aj) ~ a g+ 2-j, and 

a g+ x + ai-b 1 , j = 1, 
b g+ 2-j + a 1 -b 1 , j > 1. 

Let the homology classes be grouped symbolically into row vectors of length g: 

g L = (ai, ...,%), b = (b 1 ,...,b g ). 

Using this notation, and using (|4^|) to substitute for a g+ i, we may express the results of the 
proposition in more compact form as 

t;(o)~oM-l, tr(6) ~6'M + a(M - *M) -I, 

where £ is the row vector (£, 0, . . . , 0), and 



M 



-10 
-10 
-10 



0" 

1 

1 



(In particular, M 



-1 
-1 1 



1 
1 1 



when g = 2.) 








Remark 4.2. It is also convenient to group the holomorphic differentials of £ into column 
vectors. For example, if the 'coordinate differentials' Uj = (\ 9 ~i / /i)dA, which are holomorphic 
for I < j < g, are grouped into a column vector u, then the vector of normalized holomorphic 
differentials is given by 

to = 2mA~ 1 u, 
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where A — u ■ a is the g x g matrix of period integrals of the coordinate differentials, and we 
represent integration of a differential over a cycle by a dot. Then the Riemann matrix B is 
given by 

B = 2mAr 1 (u-b). 



Using the formula (f4*H|) and the transformation of the a-cycles under a, it is easy to check 
that the Abelian differentials transform as follows: 

a*dQi = — dfli, a*dQ2 = d^2, cx*df23 = dQ% — u\. 

(For example, a*dfls has the same asymptotics as df2 3 , but has period —2m on the cycle 
a x .) Integrating each of these over the 6-cycles, we obtain 

MV = -V, MW = W, Mr = r + (B - 27ril)i, (45) 

where the subscript one will denote the first column in a matrix or first entry in a vector. 
Similarly, a*uj = Mu, giving 

B = MB t M + 2ni(I - M f M), (46) 

Note that M 2 = I, and V, W belong respectively to the —1 and +l-eigenspaces of M, 
which we will denote as M.^ and A4 + . It is easy to verify that dim.M- = \g/2~\ while 
dimM+= \{g-l)/2\ 

Since the symmetric assumption c = Yl^=i = 0> then from the appendix gives 

1 J 1 \ CI 1 
E = > / XuOj = tr (\lu ■ a) = tr (a*(Xu) ■ a(a)) 

7T1 ^ L J 7T1 7T1 

= -1 tr (a*(Ac^) ■ (aM - £)) = - {XMU ' ~ M) + 2tr (cr*(AA -1 ^) ■ *) 
7ri 7ri 

tr(AM 2 cj-a) ^A,.^ f,f\ 9+1 - j d\ 



7T1 



(In the last step, the integral has a residue about oo + only when j = 1.) Therefore, using 

E = 2-Ki{A- 1 ) ll = \V l . (47) 

4.2. Results. The following results apply only to filaments reconstructed at the origin, i.e., 
A = in the Sym-Pohlmeyer formula. Note that, with this choice, smooth closure of the 
filament is not automatic. When g is even, the origin is automatically a zero of dOi, but it 
is not automatically a double point. 

Theorem 4.3. Assume E is symmetric. Let q(s, t) be the NLS potential associated to (E, D), 
given by (jB^j) . Let A C Z 9 be the one- dimensional lattice of vectors m such that mi = and 
m 2 = . . . = m g is odd. (Note that A C M.+-) Suppose t is such that 

D + MD - 2iWt = Trim (48) 

for some m G A. Then q(s,to) is the potential for a planar filament. 

(Note that, following ^5.21 the vector D may be assumed to be purely imaginary.) 

Corollary 4.4. If E is symmetric and g — 2 or g = 3, then M. + is one-dimensional, and 
the filament is planar at regular intervals in time. 
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By observing that, in low genus, W must be a multiple of *(0, 1, . . . , 1) G A4 + , we deduce 
that the gap between planar times is Tr/W g . It is possible that higher-genus filaments may 
be planar at isolated times, but Thm. 14.31 only implies this if the projection of D into A4+ 
lies along one of a countable number of lines parallel to W. If that is the case, then planarity 
at repeated times is possible if W is a multiple of *(0, !,...,!). 



Proof of Theorem \4-3[ Let D' = D — iWt . To prove planarity at time t it suffices to prove 
that qo(s), which we define as 

g(s,t ) giVg + r-DQ 

qo{3) '= Aexp(iivto) = eM ~ lEs) 0(iV 8 -DO ' (49) 

is the potential of a planar curve, i.e., the real and imaginary parts of go, which are the 
curvatures for a natural framing along the curve, satisfy a constant-coefficient homogeneous 
linear equation (see, e.g., |2j). To this end, compute 

0(-iVs + f + D') 0(-iVs + r + D') 



q (s) = exp(+iEs) = exp(iEs 



0(-iVs + D') ^ v ' 0(iVs-D') ' 

using the fact that D' is pure imaginary, f = r mod 27ri, and is an even function. Then, 
using D' = —MY)' + 7rim and (J45|) . we have 

— ^ ^ N 0(M(iVs + r + B l -D' + Trim)) 
qo (s) = exp(iEs) ms _ W) . 

The theta-function in the numerator is given by the series 



n£Z3 



exp(n, M(iVs + r + B x - D' + Trim) + \Bu) 

= Y exp(*Mn, iVs + r + B x - D' + Trim + \MBn) 

n 

= J2 exp(*Mn, iVs + r + B x - D' + Trim + \{B l M + 2tt1(M - *M))n) 

n 

= ^exp( i Mn, iVs + r + £>i — D' + \B 1 Mvl) exp(*Mn, Tri(m + (M — *M)n)), 



where in the third line we use the substitution MB = B l M + 2ni{M — t M), derived from 
Now we observe that the second factor in the last line is equal to one. For, let s = n 2 + . . .+n g . 
Then 

(Mn, m + (M — *M)n) = sm 2 + s(— n\ — s) — riis = s(rri2 — 2ni — s). 

When s is odd, the factor (m 2 — 2n\ — s) is even, and vice-versa. Hence, because *Mn ranges 
over all points of Z 9 , 

0(M(iVs + r + B x - D' + Trim)) = 0(iVs + r + B± - D'). 

Next, using the quasi-periodicity of 9, 

0(iVs + r + Bi - D') /• T7i \ . / . T r n/ J(iVs + r-D') 



g (s) = exp(iEs) g (iVs _ D ,) = exp(iEs) exp(-iK lS - n + D ly ^ _ ^ 

/ / ,0(iVs + r-D') . , , , N 

= exp(-iEs) exp(D 1 - n) -— r; _ — = exp(D 1 - ri)g (s), 

using the formula f|47|) for E. Since go( s ) is a constant multiple of its complex conjugate, it 
is the potential of a planar curve. □ 
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When a filament is planar at regular intervals (as predicted by Cor. 14. 4|) . these planar 
forms are congruent to each other. For, if t is a planar time and we increase time to 
t\ = to + 2ix/W g) then D' is decreased by an integer vector times 27ri. The modified potential 
go, defined in (JH2J), is unchanged. So, q(s,ti)/q(s,to) is a unit modulus constant, and the 
filament at time t 2 is congruent to that at time to- Moreover, the isometry of M 3 that carries 
the filament at time to to time t\ will be repeated in the evolution from t\ to subsequent 
planar times t 2 ,t3, etc. However, computer experiments in genus two show that the curve 
at time ti/ 2 = t + 7r/W ff is also congruent to the curve at time to, but this doesn't happen 
in genus three. 

Remark 4.5. Another phenomenon (which so far is only experimentally verified) is that 
in genus two, between the planar times to, ti/2, ti, . . ., the curve lies on a sphere whose 
center and radius vary with time. (This can be checked numerically by evaluating the real 
and imaginary parts of the potential q(s, t) at a fixed time, and verifying that they satisfy 
an inhomogeneous linear equation.) Vortex filament flow solutions obtained as Bcklund 
transformations of translating circles also lie on spheres ||. 

Proposition 4.6. Let q(s) be the potential for a curve 'y(s) in M. 3 . If q(—s) = q(s), then 
the curve is symmetric under reflection in the normal plane through 7(0). If q(s) = q(—s), 
then the curve is symmetric under a 180-degree rotation about the normal line through 7(0). 

Proof. Let T(s), Ni(s), N 2 (s) be an oriented natural frame along 7, with natural curvatures 
ki, k 2 which are the real and imaginary parts of q(s) respectively. Let 7(5) = 'y(-s). Then 
in the first case (— T(— s), AT 1 (— s), N 2 (— s)) is an oppositely-oriented natural frame for 7 
with the same curvature functions as the given oriented frame for 7. Hence 7(s) must be 
congruent to 7(s) under an orientation-reversing isometry. Since that motion must fix 7(0) 
and iVi(O), AT 2 (0), then it must be reflection in the normal plane through 7(0). In the second 
case, (— T(— s), Ni(—s), —N 2 (—s)) is a similarly-oriented natural frame for 7 with the same 
curvatures, and 7(s) must be congruent to 7(s) under an orientation-preserving isometry 
that fixes 7(0) and the line through 7(0) along iVi(O). (Since k 2 (0) = 0, this is the normal 
line at 7(0).) □ 

Proposition 4.7. Let X be symmetric, and assume there are sq and to such that 

D - (iVs + iWt ) = vrik (50) 

for some vector k of odd integers. Then the filament is symmetric under a 180-degree rotation 
about the normal line through 7(so^o) 

Proof. Let s' = s — sq, and let 

, q(s,t ) . .0(iW + r-7rik) 

Aexp{— iEsq + iNto) v{ivs' — 7rikJ 

It is easy to see that 

, „ 9(-iVs' + r + vrik) . „ 6(-iVs' + r - 7rik) , /x 

qoW = exp( +1 ^) \_ iVs , + A) > = expos') \_ iVa ,_ A) > = »(-,'), 

using the 27ri-periodicity of 9. Since |A|go(s') is a potential for the filament 7(s, to), then the 
result follows by Prop. 14.61 □ 

Theorem 4.8. Assume S is symmetric and of genus 2. Then there is a fixed plane in R 3 
such that at each time the filament is symmetric under reflection in that plane. At each 
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'planar time ' t given by Thm. \4'J[ the plane containing the filament is perpendicular to the 
plane fixed by the reflection. 

Proof. Let to be a planar time, i.e., suppose (J48j) is true for some m e A. Then there exist 
s and an odd integer vector k such that (|5U|) is satisfied. For, (f4"Sj) says that the projection 
of D' = D — iWt into Ai + is |m. Given this, then (J5Uj) implies that k + Mk = m, which 
amounts to 2k 2 — k\ — mi- Given any k that satisfies this condition, there exists an sq that 
satisfies (|5UJ). because V spans Ji4-. 

Thus, the curve is symmetric under rotation T about an axis L through "f(so, to). However, 
since the filament is also planar at time to, the rotation T can be replaced by the reflection 
1Z in the plane through L perpendicular to the plane containing the filament. 

To complete the proof, we use the symmetry properties of the vortex filament flow 

It = Is X Iss- 

If ~f(s,t) is a solution, and T is a orientation-preserving isometry of M 3 , then it is easy to 
check that 7(5, t) := T( , y(s, t)) is another solution. In particular, if 7(3, t ) = j(s — c, t ) for 
some c and some time to 5 then this equation holds at all times. In other words, symmetry 
under T is preserved by the flow. 

Similarly, if 1Z is an orientation-reversing isometry, then 7(5, t) := TZ("y(—s,t)) is another 
solution. Since 7 is symmetric under reflection 1Z at time to, and this symmetry reverses 
orientation along 7, then TZ("y(—s,t)) = 7(3 — c, to) for some c, and it follows that this 
relationship persists in time. □ 

Remark 4.9. At a planar time, a filament must cross the axis of symmetry at least once, 
since the filament is connected to its reflection. If it crosses the axis in a direction that is not 
perpendicular to the axis, a self- intersection occurs at this point. The fact that reflection 
symmetry persists as the curve lifts out of the plane implies that these self-intersections 
also persist (see Figures and EJ) . However, it is possible that, as the curve evolves, two 
self-intersections may coalesce and disappear, while maintaining reflection symmetry. 

Note that not all planar curves that are symmetric about an axis have self-intersections 
which are preserved by the flow. It is easy to draw examples where the reflection symmetry 
preserves the orientation along the curve (rather than reversing it, which is the property 
used in the proof of Theorem 14. 8|) . Then, under the filament flow, the strands immediately 
pull away in opposite directions from the self-intersections. In this case, the symmetry is 
the restriction of an orientation-preserving 180-degree rotation in M 3 , and such rotational 
symmetry persists as the filament evolves. 

Remark 4.10. The argument in the proof of Theorem 2 relies on the fact that, when g = 2, 
V and W span M. 9 . However, given a symmetric hyperelliptic curve of genus g > 3, there are 
plenty of potentials to which the argument in the proof of Theorem 14.81 applies (see Figure 
EJ). For, given any vector k of odd integers, we may choose an imaginary vector D so that 
(J5U|) is satisfied for some s and t . Then D also satisfies (J4"S|) if and only if the differences 
k g _j — k 2+ j, j > 0, are all equal. The latter condition is vacuous when g = 3. Thus, in genus 
3 the set of appropriate divisors corresponds to the image of the affine lattice of odd integer 
vectors in the quotient vector space M 3 /{V, W}. (Changing D by a multiple of V or W 
merely corresponds to a shift in so or to-) 



28 



A. CALINI AND T. IVEY 






Figure 7. Genus 2 symmetric filament, generated using D = 0, at times 
t = (Ak — l)(7r/2W / 2), for k from to 5. This illustrates the periodic planarity 
of Thm. 14.31 the sphericity of Remark |4.51 and persistent planes of symmetry 
of Thm. EH 
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Figure 8. Genus 3 symmetric filament, generated using D = 0, at times 
t = {Ak — l)(7r/2W / 2), for k from to 5.. Note the persistent self-intersections 
in the plane of reflection symmetry. The curve returns to its original shape, 
in a rotated plane. 
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5. Appendix 

In this appendix, we will indicate how the formulas for finite-gap NLS Baker-Akhiezer 
functions developed in ^j21are related to those given in the monograph by Belokolos, Bobenko, 
Enol'skii, Its and Matveev pQ. We will also develop (from a slightly different perspective 
than pQ) the reality condition which must be satisfied by the divisor, and discuss how to 
compute some of the quantities used in the expression for the finite-gap NLS potential. 

5.1. Alternative Eigenfunction Formulas. From §2.1\ we had 
^(P; s,t) = ^^exp (i S (fi!(P) - f ) + it(fi 2 (P) + f )) 



6{A{P) + iVs + iWt - A(V+) - K) 9(A(oo_) - A(V+) - K) 
6(A(oo-) + iVs + iWt - Ap+) - K) 6(A(P) - A{V + ) - K) 



and 



^ 2 (P; s, t) = ±¥L (ia(ni(P) + f ) + i*(fi 2 (P) - f )) 

6{A{P) + iVs + iWt - Alp J) - K) 6(A(oo + ) - Aip_) - K) 
X 6{A{oo + ) + iVs + iWt - Alp J) - K) 6(A(P) - Alp-) - K) ' 
Recall that is the Abelian differential on S, with zero a-periods, such that 

dfi 3 ~ ±— - as A — * oo±. 
A 

Then ^df^ = 27ri, where I is the clockwise cycle around oo + shown in Figure ^ 

Lemma 5.1. § h dVt^ = —Tj, where Tj denotes the jth entry of vector r. 

This will be proved at the end of ^5.31 For now, we note that the corresponding Abelian 
integral Q3 behaves like ±log(A) near oo±. Furthermore, we observe that the quotient 

pv sv e(A(P)-D) 

is a well-defined meromorphic function on S. For, if the paths of integration in the integrals 
f2 3 (P) and A(P) are changed by the addition of the same a-cycle, each factor in /(P) is 
unchanged; if the cycle bj is added, then f(P) is modified by the factor 



exp dfl-^j exp(r 3 ) = 1. 



Note that if, say, we confine the path of integration for Q 3 to the interior of the cut 
Riemann surface E described in £ 12.11 then exp(f2 3 (P)) has a pole at oo + and a zero at oo_. 
Consequently, /(P) has pole divisor T> + + oo + and zero divisor T>- + oo_. It follows that 

exp(n,(P)) \ { A ( p) _ B) J x^| (51) 

is a nonzero constant. 

We can express this constant in terms of the asymptotic behavior of Q3 and g + . Namely, 
suppose P is near oo + and lies over the real axis in the A-plane. Let Cp be a path from P 
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to l(P) in the interior of S , such that Cp is invariant under sheet interchange. Let a g+1 be 
the cycle around the rightmost branch cut, which satisfies the ho mo logical equivalence (J44j) 
within £ — {oo + , oo_}. Then we have 



dfi 3 = / dn 3 = dtt 3 = -2vri + / dQ 3 , 

Cp J r(Cp) J Cp-\-a,g-\-\ J Cp 

where r : (A,/i) \- > (A, fx) is complex conjugation on E. Therefore, 

lim (/ dft 3 -21ogA(P) ] =7ri-log/5 (52) 

P->oo+ \J Cp J 

for some positive real number f3. (In terms of the notation in £Q, the above limit is written 
as — logc^o for a negative real number ujq.) Therefore, 

exp(fi 3 (P)) ~ -J=A as P -> oo + . 

Suppose that g+(P) = — - + 0(A~ 2 ) when P is near oo + , for some constant a. Then letting 

aX 

P tend to oo + in (JoTj) gives 

exp( n3 (/-)) g( ^ P l; D ^ r) x 9+(P> - 1 " <D) 



S{A(P)-T>) gJP) a% fpe(D-rY 
Solving this equation for g~(P) and substituting in ip 2 gives 

^ = g+ (P)^exp (is(n a (P) + f ) + i*(0 2 (P) - f ) + Os(P)) 

0(.4.(P) + iVs + iWt — D — r) 0(D 

X 



6(i~Vs + iWt - D) 6(A(P) - D) 

which agrees with the formula 4.1.16 for ip 2 in 1J, up to the multiplicative factor of g + (P). 
Similarly, if) 1 coincides with the formula for ip 1 in pQ except for multiplication by g + (P). 
Since our vector-valued function and that of pQ differ only by multiplication by a scalar that 
is independent of s and t, both satisfy the NLS linear system for the same potential q(s,t). 

5.2. Reality Conditions. By inserting the matrix \1/ defined by (fTTj) into the NLS linear 
system, and expanding near oo_ in powers of A, we obtain 

, , , , -rn . Ar N #(iVs + iWt - D + r) , A 26(B) 

q (s,t)=AeM-iEs + m { ms + mt _ D) > where A = -jjl-L, (53) 

and 

r 4/3 , 9(i\Fs + iWt — D — r) , . 

9M) = ^exp(i gs -iiW) V^iWt-D) ■ (54 > 

It is not immediate that the second expression is the complex conjugate of the first; in 
fact, the reality conditions on D come from imposing this condition. Before computing this, 
we need to calculate the imaginary parts of the members of ([53)1 . Since T*dQj = dQj and 

r{h) ~b k + e + Y, a i 

within £ — {oo + , oo_}, then V and W, being the 6-periods of dOi and df^, are real vectors. 
Similarly, Lemma 15. II implies that 

rk = n + 27ri. 
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Moreover, since r*u; fc = u7/7, then the Riemann matrix B, defined by satisfies 

( for k = I 

Bkl = \B kl + 2vri for k ^ I. (55) 



It follows that 9(z) = 9(z). 

Now comparing (J53j) and (p^fj) gives the reality condition 

-j0(-iVs - iWt - D + r) 4/3 #(iVs + iWt — D — r) 



0(-iVs - iWt - D A 0(iVs + iWt-D) 
Since poles and zeros must match up, we must have 

D + D = 27rin + Bm, (56) 

where m, nGZ 9 are such that 

|A| 2 exp((m,r)) = 4/3. (57) 
Taking the imaginary part of (|56|h with (jB3|) taken into account, gives 

2nj = ^^mfc. (58) 

(In particular, n = when g = 1.) Thus, all the entries of m have the same parity, even or 
odd. When g is even, summing (|58|) over j shows that m must be all even. If g is odd and 
m is all odd, the fact that \m(rj) = it would imply that the left-hand side of (JH7|) is negative. 
Thus, we can assume that m is all even. In this case, one can let D' = D — ^Bm; then D' 
is purely imaginary, and (}53|) . (J5?]) are still valid with D replaced by D' and m replaced by 
zero. We now do this, and then our reality condition amounts to the vector D being purely 
imaginary: 

D = -D. (59) 

We can, without loss of generality, assume that A is real and positive; then (|H7jl gives 
A = 2y/]3. From (J5H|) we then get ot.y/~f3 = 29(D)/ 9 (D — r). So, omitting the common scalar 
multiples g + (P)9(D) / 9(A.(P) — D) from the eigenfunctions gives the simplified formulas 

9(A(P) + iVs + iWt - D) 



^ = exp {is(Q l (P) - f ) + it(Q 2 (P) + f )) 
^ = -iexp (is^P) + f ) + it(Q 2 (P) - f )) 



9(iVs + iWt - D) 
9(A(P) + iVs + iWt - D 
9(i~Vs + iWt - D) 



We will now interpret the reality condition (J5U|) in terms of the divisors on E. 

Proposition 5.2. IfD = -D, then V_ = lt(D + ). 

(Recall that i denotes the sheet interchange automorphism and r denotes complex conjuga- 
tion on E.) 

Proof. The function f(P) = 9(A{P) — D) has zero divisor T> + . (To make / well-defined, we 
restrict to the cut Riemann surface Eo, assuming that the homology basis avoids the points 
of T>+ and T>_.) Therefore, the function 



f(LT(P)) = 9[A(ir(P))+D) 
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has zero divisor lt(V + ). But since oo + = ir(oo_) and (tr)*ujj = —uJ], then 

A(lt(P))= / LU+ / LU = r-A(P). 

Joo- J lt(oo-) 

Using the fact that rj = rj modulo 27ri, 

6 (A{it(P)) + D) = 9(f - A{P) + D) = 9{A{P) — D — r) = d(A{P) - A{VJ) - 7C), 
which has D- as its zero divisor. □ 

Proposition 5.3. Assuming T> satisfies the Genericity Condition o/3B then T>_ = lt{T> + ) 
if and only if lt(T>) ~ T>, where ~ denotes linearly equivalence of divisors. 

Before proving the lemma, we recall some pertinent facts about divisors on a Riemann 
surface (see, e.g., ^2] or [221 )• For any divisor T>, the meromorphic functions / such that 
T> + (/) > form a vector space C(T>) with dimension l(T>). For example, a positive divisor 
always has l(V) > 1, because of the constant functions. If two divisors V,V are linearly 
equivalent, then l(V) = l(T>'). 

The Genericity Condition can be stated as 

l(V — oo_ — oo + ) = 0. (60) 

This implies that 1{T> + ) = l(T> — oo + ) = 1, and similarly for X>_ and T> — oo_. For, if 
l(T> — oo + ) > 1, let fx, f 2 G C(V — oo + ) be linearly independent; then there exists a nontrivial 
linear combination cif\ + C2/2 that vanishes at oo_, and this contradicts (J6Uj) . 

Proof of Prop. \5.c\ Suppose that T>_ = lt(T> + ). By construction, 

V + = V-oo + + (f) 

for some meromorphic function /. Appying it to both sides gives 



O IT) . 

Since ZL ~ X> — oo_, then X> ~ lt(T>). 

The same computation, run backwards, shows that if T> ~ lt(T>) then £L ~ 6r(X> + ), i.e., 

D_ + (/) = ir(X> + ) 

for some meromorphic function /. However, l(T>_) = 1 implies that T>_ + (/) > only when 
/ is a constant. Thus, XL = lt(T> + ). □ 

5.3. Computing E, N, r and (3. Suppose functions Q, Q' are holomorphic on the cut 
Riemann surface S . Then by Stokes' Theorem, 

0= / d(Q!dn)=<f n'dtt=Y^A / j B j - B'jAj, (61) 

JT, JdT,o j = i 

where Aj = § Q, Bj = § b .Q, and similarly for the primes. (See, for example, ^2]-) We 
will use this relationship to derive formulas for E and iV and to prove Lemma 15.11 (Our 
formula for E is equivalent to one given without proof in [I].) The differentials we use will 
actually have singularities at oo±, so the left-hand side of (JBTj) will be modified by adding 
the integrals of Q'dfl along clockwise cycles about these points: 

= V A'jBj - BjAj + I Q'dQ. (62) 
j= i Jt+^W 
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Each of these extra terms will be — 2ni times the residue of the differential at the singularity. 

First, let Q' = Qi and dQ = (X 9 /fi)dX. By expanding fi 2 in terms of A and applying the 
binomial series, we get that, near oo±, 

X 9 + 

where c = Y^f=i <\j an d 



±(l-|A- 1 -rfA- 2 + 0(A- 3 )) 



2<y+2 

c 2 "2£a?) (63) 



This, together with 



Oi = ±(A+| + o(l)) 



allows us to expand the integrand near oo±: 



van = -A 2 ( i - ( ^y^) a- 1 + o(A- 2 ) ) d(A- 1 ). 



Since dTl' = dOi has no a-periods, from (Jozjl we have 

9 



f X 9 

= -2ni(E - c) - V o f — dA - 



Thus, 

Next, let Q' = f2 2 instead. Then 

n'dn = -A 3 ^2 + cX- 1 + (y + 2rf + + °( A ~ 3 )) d ( A_1 ) 

near oo±, and (|U2*j) gives 

= -2vri (JV + Ad + c 2 ) - ^Wjl A 



so that 



-dA, 
/i 



— YwJ -dA. 



AT = -4rf - c 2 — dA. (65) 



We note that, by letting df2 = and Q' = Q± or Q 2 , similar calculations give 

Q = 2c jl -V j , = Ac j2 + 2cc jl -W j , (66) 
where the constants such that 

9 \g-k 



y^Cjk dA. 



fc=l 



Thus, computing just the a-periods of the coordinate differentials u k = dA for < k < g, 

we can determine E, N and the 6-periods of dfti, dQ 2 algebraically. 
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In most experiments we have carried out, it is straightforward to calculate the a-periods of 
the coordinate differentials by numerical integration. Numerical techniques are also practical 
for calculating (3 in the limit ()52j1 . For this purpose, we can rewrite that limit as 

27ri-log/? = 21og(A 23+2 ) + 2 lim / (dfi 3 - A~ x dA) , 

P^oc+ J Tp 

where V p is the portion of Cp from A 29 +2 to P on the upper sheet. Now the improper integral 
converges. 

Proof of Lemma \5. 11 If we let Q' = Q% and dQ = Uj, we must make additional cuts. Let 

Eq be the cut Riemann surface So (shown in Figure El) with an additional cut from oo_ to 

oo + , along II U t(n) . Then ^3 is well-defined and holomorphic on E , with its values of ^3 

differing by 27ri on each side of this cut. From (JoT|) we obtain 

r r°°+ 
= -27ri <f> dtt 3 - 2vri / ujj. 

J bj J oo_ 

□ 
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